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We investigate the transport of excitations through a chain of atoms with non-local dissipation 
introduced through coupling to additional short-lived states. The system is described by an effective 
spin-1/2 model where the ratio of the exchange interaction strength to the reservoir coupling strength 
determines the type of transport, including coherent exciton motion, incoherent hopping and a 
regime in which an emergent length scale leads to a preferred hopping distance far beyond nearest 
neighbors. For multiple impurities, the dissipation gives rise to strong nearest-neighbor correlations 
and entanglement. These results highlight the importance of non-trivial dissipation, correlations and 
many-body effects in recent experiments on the dipole-mediated transport of Rydberg excitations. 


The transport of energy, charge, or spin is of funda¬ 
mental importance in diverse settings, ranging from the 
operation of nanoelectronic and spintronic devices m, 
to the dynamics of electron-hole pairs in organic semi¬ 
conductors [3IE] and in natural processes such as pho¬ 
tosynthesis [ 31 ]. In these systems different elementary 
excitations and basic transport mechanisms can give rise 
to very different behaviour including coherent exciton 
motion, thermally activated diffusion or even collective 
fluid-like dynamics [TO]. However, understanding or ex¬ 
ploiting these differences (e.g. for applications in pho¬ 
tovoltaic devices) is extremely challenging, as they de¬ 
pend on the complex interplay between quantum statis¬ 
tics, coherence, confinement, disorder, and the nature 
of the interactions between the constituent particles. A 
key question is how dissipation, through the coupling to 
reservoirs, leads to a cross-over between coherent and in¬ 
coherent motion. While dissipation is usually assumed to 
destroy coherence, it is becoming evident that certain dis¬ 
sipative processes, noise or specially structured environ¬ 
ments may even enhance coherence or enable new modes 
of transport, both in engineered quantum systems and in 
natural ones [Hiiniin]. 

Here we analyse the exciton-like motion of individual 
excitations through a network of dipole interacting quan¬ 
tum systems (i.e. atoms, molecules or quantum dots) in 
the presence of a specially engineered reservoir (Fig. Hb- 
Each subsystem is assumed to be coupled to one or more 
additional short-lived states. The populations of these 
states are determined by the coupling Adds and by their 
proximity to any excitations, which provides a handle to 
introduce and control new types of dissipation in the sys¬ 
tem. By eliminating the short-lived states we show that 
the resulting system is characterised by coherent hopping 
and non-local dissipative terms which lead to correlated 
exciton-mot ion. Simulations of the single excitation and 
the two excitation dynamics demonstrate this has a dra¬ 



FIG. 1. (Color online) Transport of excitations (red spheres) 
through a chain of optically-dressed-atoms with site index j 
(green spheres). (a) Excitonic motion arises through coherent 
long range exchange interactions J and through dissipative 
processes mediated by a tailored reservoir coupling, (b) Each 
site is represented by an atom which can be in either the im¬ 
purity Rydberg state \p) or the auxiliary states \g), \e) or |s). 
Dipolar interactions result in exchange between \s) and \p) 
Rydberg states. States \g) |e) and |e) |s) are coupled 

with Rabi frequencies k and D respectively. The state |e) 
spontaneously decays to the \g) state with the rate Fe. The 
state space He comprising |e) or |s) excitations can be adia- 
batically eliminated resulting in an effective spin-1/2 model 
for the impurity and dressed states. 


matic effect on the transport properties, leading for ex¬ 
ample to the emergence of a new length scale for hopping 
and strong spatial correlations even at steady state. 

As a physically realizable system we propose a chain 
of ultracold atoms in which Rydberg states with large 
transition electric dipole moments are optically coupled 
to low-lying electronic states, which provides a clean 
synthetic system to investigate energy transfer dynam¬ 
ics and transport including controllable interactions and 
dissipation Recently the exciton-like migra¬ 

tion of an ensemble of Rydberg impurities in a bulk 
system was observed using electromagnetically-induced- 
transparency (EIT) of a background atomic gas as an am¬ 
plifier [25]. Remarkably, the flgure-of-merit for exciton 
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migration, the diffusion length Dr ^ 50 — 100 jam 

(with D the diffusion coefficient and r the typical lifetime 
of the Rydberg state) was an order of magnitude larger 
than in even the purest organic semiconductors [26] , and 
can be tuned by the probing light fields which act as 
a controllable environment. In another experiment the 
coherent hopping of a single excitation amongst three 
equidistantly spaced Rydberg atoms was observed m 
Future experiments will be able to probe the coupled co¬ 
herent and incoherent motion with high spatial and tem¬ 
poral resolution, while the strength and nature of the 
dipole coupling, the degree of disorder, excitation den¬ 
sity and the role of the environment can all be controlled, 
allowing unprecedented opportunities to investigate the 
fundamental processes at play. 

Consider a system of ultracold atoms initially prepared 
in either a Rydberg state with orbital angular momentum 
I = 1 (impurity state \p)) or the electronic ground state 
1^). The 1^) state is weakly optically coupled via an FIT 
resonance to a short-lived excited state \e) and an / = 0 
Rydberg state \s) (Fig. The 1^) ^ |e) probe tran¬ 
sition is driven with Rabi frequency k and the \e) js) 
coupling transition is driven with Rabi frequency Q with 
hz ft. The state |e) spontaneously decays to \g) with a 
rate Fg, while the other states are assumed to be stable on 
the time scale of the dynamics. Migration of the \p) exci¬ 
tations among the dressed ground state atoms occurs via 
the \s) state admixture (by \sj)\pk) ^ \Pj)\sk) exchange, 
with j and k the site indices) [ 20 ], while the |e) state ad¬ 
mixture introduces a controllable environment [28] , both 
of which are influenced by the competition between the 
EIT laser fields and the Rydberg-Rydberg interactions. 

The Hamiltonian describing the system is given by 

{h = i) 

H = ( 1 ) 

where Hq = \^j^\sj){ej\ + + h.c. ac¬ 

counts for the single-atom laser couplings and 
= Cm/\xj — Xk\^. Here we mainly consider 
dipolar interactions corresponding to m = 3, however 
m = 6 can also be realised using nonresonant van der 
Waals interactions between Rydberg states [29H3T]. For 
simplicity the interaction coefficient Cm is assumed to 
be independent of j and k. Additionally, spontaneous 
decay of the intermediate excited states is included 
through Lindblad terms leading to the master equation 
for the density matrix p 

p = -i[H,p\ + C[p\ (2) 

where C[p] = LjpVj — ^ [L^-Ljp-\- pV-Lj) and each of 
the Lindblad operators represents a single decay channel 
{Lj = \/Tl\gj){ej\). 

Simulating the open-system dynamics for more than 
approximately five four-level atoms is beyond what can 


be readily performed using exact numerical methods (e.g. 
Monte Carlo wave-function techniques). However, in the 
weak probe regime most relevant to experiments {n 
11 , Fg) the populations of the |e) and \s) states are always 
small. Furthermore the dynamics between these states 
due to the laser coupling is fast compared to the hopping 
dynamics. Therefore the many-body states which include 
|e) or \s) can be adiabatically eliminated [32f[34] . 

For this we use the effective operator approach [35] . 
which allows for a simple interpretation of the dominant 
processes in terms of coherent and incoherent coupling 
rates. First the state space is separated into a slowly 
evolving ‘ground-state’ subspace involving only \g) , \p) 
states and a rapidly evolving ‘excited-state’ subspace in¬ 
cluding |e), \s) states (Fig. &)• The ground state Hamil¬ 
tonian contains no direct couplings Hg = 0 , while the 
coupling laser and the dipole-dipole interactions enter the 
Hamiltonian for the excited state manifold ^g, and the 
probe laser weakly couples ground and excited subspaces 
through Vj^ = Vl. The master equation governing the 
evolution of the ground states is then defined through the 
operators = —'F_Re[x]+^^ and where x 

is the solution to the matrix equation Hnh^ = with 
the non-Hermitian operator Hnh = ~ (V^) 

To simulate the time evolution of the system for a given 
set of parameters we first numerically obtain x in order to 
define an effective master equation [as in Eq. (© but us¬ 
ing the operators and which can then be solved 
in the usual fashion. 

Single exciton dynamics:- In the special case of a sin¬ 
gle \p) excitation and by neglecting beyond second-order 
interactions (corresponding to couplings between states 
involving more than one dressed-atom), simple analytic 
expressions for the effective operators can be found. The 
resulting A^-site effective master equation can be written 

+ h.c. 

k>j 

Lf = Y.^^fCdC)SXS'l. - Cidhk)SXsi (3) 


where dj^k = \xj — Xk\ and we have made use of the 
spin raising and lowering operators (5+ and S- respec¬ 
tively with = \pk) {gk\ where \g) is the dressed ground 
state). 

The resulting effective operators involve three 
terms: (I) effective coherent exchange interactions 

incoherent hopping T = 

and (3) irreversible dephasing 7 = acting 

on each site. Here K = n^/Te and V = (Rc/dj^k)^ 
where Rg = ( 2 rgC 3 /lI^)^/^ is the dipole blockade ra¬ 
dius [25] [36] . These two parameters have simple interpre¬ 
tations as the two-level-atom photon scattering rate, and 
the dipole-dipole interaction energy scaled by the EIT 
bandwidth respectively. K~^ defines a natural timescale 
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for the dynamics whereas V{dj^k) is responsible for the 
competition between coherent exchange, hopping and de¬ 
phasing. Equations can be used to efficiently simu¬ 
late the open-quantum-system dynamics of a single impu¬ 
rity immersed in a background of hundreds of optically- 
dressed atoms in arbitrary geometries. It is interest¬ 
ing to note that this effective model is closely related 
to the widely used Haken-Strobl-Reineker (HSR) model 
for exciton motion in the presence of noise EH EH], in¬ 
cluding two distinct decoherence mechanisms which orig¬ 
inate from the spontaneous decay of the |e) states. How¬ 
ever, the incoherent hopping jump operators have an un¬ 
usual form: describing correlated jumps 

from states \pj) to collective states involving many neigh¬ 
bouring sites k [34]. This is in contrast to the original 
HSR model which ignores the influence of the system 
on the reservoir and assumes uncorrelated fluctuations 
e.g. Therefore, the system studied 

here offers the possibility to use dissipation controlled 
via EIT resonances to control both the strength and na¬ 
ture of decoherence in these systems which, for example, 
is an active ingredient of proposals for quantum state 
preparation by reservoir engineering [391143] . 

Fig-E a) shows the rates for coherent exchange, hop¬ 
ping and dephasing terms as a function of the relative 
distance between two sites. Eor short distances coher¬ 
ent exchange is suppressed while at large distances it is 
determined by the dj^ dependence of the dipolar interac¬ 
tion J « /VP‘)Vdd‘ The peak exchange rate occurs for 
dj^k = Re where J = K/A. At this distance however in¬ 
coherent hopping and irreversible dephasing are equally 
important. Eor dj^k 0 dephasing dominates, saturat¬ 
ing at a rate given by the two-level atom scattering rate 
-f = K. 

We now turn to the analysis of spin dynamics in 
this system for the case of atoms arranged in a one¬ 
dimensional chain with intersite separation a (as can be 
produced for example in an optical lattice). We antic¬ 
ipate three main regimes: (1) for V{a) 1 coherent 
exchange dominates over decoherence, (2) for V{a) ^ 1 
decoherence becomes important leading to classical hop¬ 
ping, (3) for V{a) ^ 1 strong dephasing suppresses co¬ 
herent exchange, leaving incoherent hopping with a char¬ 
acteristic hopping distance comparable to Rc. Simula¬ 
tions for a single impurity for A" = 121 sites are shown in 
Eigs. I^b-e). The strength of the dipolar interactions is 
varied through the nearest neighbour Rydberg-Rydberg 
interaction coefficient V{a) = R^/a^. We calculate the 
impurity probability distribution Pj{t) = TT{njp{xjR)) 
(where fij = Si^S^_) on each site for different values of 
V{a). 

Eor V{a) <C 1 the dynamics are characteristic of a 
quantum random walk, with light-cone-like spreading 
and interference fringes in the spatial density distribution 
(Fig. 1^). In this regime coherent dipolar exchange dom¬ 



Rc) 




FIG. 2 . (Color online) (a) Distance dependence of the ef¬ 
fective coupling rates J (solid line), F (dashed line) and 7 
(dotted line) for a single impurity with m = 3 in units of the 
two-level atom scattering rate K = tP/Te. (b) Mean square 
displacement of the propagating excitation (x^) for a chain of 
A = 121 sites obtained from simulations for three values of 
the dimensionless interaction strength: V{a) = R^jP = 50 
(triangles), V{a) = 1 (circles) and V{a) = 0.02 (squares). 
The solid lines show the approximate analytic scalings in the 
three regimes described in the text, (c-e) Zoomed in density 
plots showing the impurity probability distribution at short 
times for (c) coherent transport with V{a) = 0.02 (d) diffu¬ 
sive transport with E(a) = 1 and (e) transport in the block¬ 
ade regime V{a) = 50. Bright colors indicate high probability 
density. 


inates leading to ballistic expansion of the wavepacket. 
For short times t < F(a)“^ the mean square displacement 
(including beyond-nearest neighbor exchange) evolves ac¬ 
cording to (x^)/a^ « (7r^/180)E(a)^A^t^ [Fig.j^b)] [44] . 
For intermediate interactions {V{a) ~ 1) the coupled co¬ 
herent and incoherent motion makes the dynamics more 
difficult to describe, however, from our simulations we 
find normal diffusion (Fig.[^) with (x^)/a^ = 2Dt where 
the diffusion coefficient D « [0.262 -f- 0.123E(a)]A is 
found by expanding around V{a) ~ 1. For V{a) ^ 1 
qualitatively new behavior is observed. In this regime 
there is a new preferred hopping distance given by the 
Rydberg blockade radius Rc ^ a (Fig. §). This can 
be understood as the presence of an impurity shifts 
the |s) states of nearby atoms which suppresses energy 
transfer. In this regime the underlying lattice geome¬ 
try becomes less important, leading for example to a re¬ 
duced infiuence of possible disorder in the atomic posi¬ 
tions. An expression for the mean square displacement 
can be found by neglecting coherent exchange and in¬ 
tegrating the incoherent hopping rate for all other sites 
(x^)/a^ = {7r/6)V{a)Kt. This scaling differs from the 
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FIG. 3. Two exciton dynamics and density-density correla¬ 
tions for r=l,Q = 2, = 0.05, m — ?> and vary¬ 

ing V{a). (a) Probability density as a function of time for 
V{a) = 1 showing normal diffusion, (b) Equal time correla¬ 
tion function g^‘^\0,k) for the same parameters. The colors 
correspond to g^‘^^ > 1 (dark red), g^‘^^ < 1 (blue) and g^‘^^ = 1 
(white), (c) Equal time density-density correlation function 
at steady-state as a function of V{a) [same colorscale as (b)]. 
(d) Nearest neighbour correlation function ^^^^(0,1) versus 
V{a) showing bunching around V{a) ~ 1. 


simple picture reported in [25], where continuous obser¬ 
vation of the system via EIT was assumed to lead to a 
Zeno-like slowdown of the dynamics. Instead we demon¬ 
strate that the system exhibits a more complex form of 
dissipation, which allows for rapid transport even in the 
strongly dissipative limit. 

Exciton-exciton correlations:- For times which are long 
compared to the inverse dissipation rate, each of the dif¬ 
ferent regimes described in Fig. exhibit diffusive be¬ 
havior and for a finite system the excitation becomes 
uniformly distributed over the chain (Fig. IT)- However, 
even at steady-state (p = 0) the nature of the trans¬ 
port may be revealed through the analysis of higher- 
order statistical properties. To demonstrate this, and 
to point out the importance of many-body effects, we 
present simulations of the effective master equation for 
two excitations, which could be experimentally prepared, 
for example, by tuning the excitation laser frequency to 
match the van der Waals interaction energy between a 
pair of p-excitations at a well defined distance [451146]. 
Time-dependent simulations for N = 16 sites for inter¬ 
mediate interaction strength are presented in Fig. ia,b) 
and solutions to the steady-state effective master equa¬ 
tion for a chain of N = 20 sites with varying in¬ 
teraction strength (with periodic boundary conditions) 
are shown in Fig. |^c,d). While the density distribu¬ 


tion evolves similarly to the case of a single excitation 
[Fig-i a)], we observe strong density-density correlations 
= c{hjhk)/{{nj){nk)) where the coefficient 
c = {2N — 2)/N normalizes to the case of precisely two 
excitations distributed over the chain in an uncorrelated 
manner [Fig. |^b)]. These correlations persist even at 
steady-state, long after the memory of the initial state is 
lost [Fig.j^c)]. Fig.j^d) shows how the nearest-neighbor 
correlations depend on H(a), exhibiting strong bunch¬ 
ing for V{a) ~ 1. For both strong and weak interac¬ 
tions the correlations are suppressed, indicating that they 
arise as a consequence of the competition between dipole- 
mediated exchange and distance dependent dissipation 
provided by the optically-dressed atoms. This dissipation 
is minimized for neighboring pairs and can be thought of 
as inducing an effective attraction between the impuri¬ 
ties, similar to effects in [miiH]. These correlations van¬ 
ish if the correlated jump operators of the effective master 
equation are replaced by operators which couple localized 
states alone. In the case of short-range interactions (i.e. 
van der Waals with m > 6) these correlations become 
even more pronounced, and for nearest-neighbor interac¬ 
tions the steady-state corresponds to the entangled pure 
state \D) = {N - \h ■ ■ -PjPj+i ■ • 

The described impurity plus optically-dressed atom 
system exhibits rapid spin transport controlled by a sin¬ 
gle parameter V{a) = 2FeC3/(r^^a^). In addition to the 
coherent and incoherent hopping regimes recently inves¬ 
tigated in a related disordered system [28], we find a 
blockade dominated regime in which a most-likely hop¬ 
ping distance emerges given by the dipole blockade ra¬ 
dius. In the limit of low impurity densities we derive an 
effective master equation that can be directly applied to 
current experiments on Rydberg energy transport in op¬ 
tically driven systems, helping to elucidate the interplay 
between coherent spin exchange and decoherence due to 
spontaneous decay of the dressed states. Similar physics 
could arise in other systems where correlated noise [49] 
or exciton-vibrational coupling [5QH53] significantly af¬ 
fects transport properties. For multiple excitations we 
observe strong density-density correlations which persist 
at steady-state indicating effective interactions between 
excitations mediated by the non-trivial dissipation. This 
suggests a novel and physically realizable experimental 
route towards dissipative entanglement creation [394143] 
and realizing exotic pairing mechanisms or exotic quan¬ 
tum phases in the many-body regime [54] [55] . 
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